Robust Mesoscopic Superposition of Strongly Correlated Ultracold Atoms 
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We propose a scheme to create coherent superpositions of annular flow of strongly-interacting 
bosonic atoms in a ID ring trap. The non-rotating ground state is coupled to a vortex state with 
mesoscopic angular momentum by means of a narrow potential barrier and an applied phase that 
originates from either rotation or a synthetic magnetic field. We show that superposition states 
in the Tonks-Girardeau regime are robust against single-particle loss due to the effects of strong 
correlations. The coupling between the mesoscopically distinct states scales much more favorably 
with particle number than in schemes relying on weak interactions, thus making particle numbers 
of hundreds or thousands feasible. Coherent oscillations induced by time variation of parameters 
may serve as a 'smoking gun' signature for detecting superposition states. 
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I. INTRODUCTION 



Quantum superpositions of macroscopically distinct 
states are important for our understanding of quantum 
mechanics [T] and carry great promise for enhanced preci- 
sion measurement techniques [2j. A conceptually simple 
example is the superposition \N, 0) + |0, N), where all N 
particles occupy either one or the other of two accessi- 
ble modes (e.g. spin orientations). Due to their inherent 
fragility, such maximally entangled NOON states [3] en- 
gineered in optics and spin systems have been limited to 
10 particles [¥]. The related mesoscopic superpositions 
of flux states consisting of 10 9 Cooper pairs observed 
in superconducting rings [5j [6] have proven more robust 
but their microscopic nature is debated [7HTU] . Proposals 
to create mesoscopic superpositions of ultra-cold atoms 
have so far entirely focused on NOON and closely re- 
lated states [TTUT4] and have not yet been realized. Such 
proposals suffer from severe limitations due to decoher- 
ence |15j and the unfavorable scaling of precision and 
time scales needed to produce these states [HI [H] • 



In this paper we propose a simple and experimentally 
accessible scheme for producing large, robust quantum 
superposition states of ultra-cold atoms. NOON states 
are so fragile because the loss of a single constituent par- 
ticle, if discriminant between either of the two available 
modes, destroys the superposition. We show that many- 
particle superposition states can be made more robust 
by making use of interactions between the atoms. The 
intuition is that correlations due to particle interactions 
spread single-particle observables over many modes. This 
disguises the origin of any lost particle and allows the su- 
perposition to survive. In addition, strong interactions 
remove degeneracies and therefore decrease sensitivity to 
environmental fluctuations. 
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FIG. 1. (Color online) Strongly-correlated atoms trapped in 
a narrow ring with a rotating barrier. The inset shows the 
energy levels of the super-current states with total angular 
momentum K — and K = Nh, respectively, as a function of 
the rotational phase Q. (dashed lines), and the lowest energy 
levels of N = 99 atoms in the Tonks-Girardeau regime with 
barrier strength b/L = 0.008-Eo (full lines). 



II. THE SYSTEM 

We apply this idea to bosonic atoms confined to a 
thin ring-shaped trap [TSl Qj|] that is intersected by a 
focused blue-detuned laser beam, which creates a poten- 
tial barrier for the atoms as illustrated in Fig. [l] One- 
dimensional Bose gases with variable repulsive interac- 
tions have already been realized in linear traps where 
the interactions were tuned by means of a Feshbach res- 
onance or by adjusting the trap geometry [201 IS]- We 
model the system of N atoms of mass M in a loop of 
circumference L by the one-dimensional Hamiltonian 
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where x = 6L/2ir is the atom's position on the circum- 
ference of the loop and g is the effective one-dimensional 
interaction strength between the atoms [52]. The small- 
est nonzero kinetic energy of a single particle 
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(2) 



provides a natural unit of energy for this system. The 
narrow barrier with strength b is rotated with constant 
tangential velocity v = Ml /(ML) along the circumfer- 
ence of the ring and the Hamiltonian (fTl) is formulated 
in the co-rotating frame of reference. Alternatively, the 
rotational phase fl can be applied in a non-rotating sys- 
tem by means of a synthetic magnetic field as recently 
demonstrated at NIST [23] . 

In the absence of the barrier (6 = 0), the angular mo- 
mentum in the ground state is quantized to integer mul- 
tiples of Nh 24. 25 . In this case, the Hamiltonian (fTl) 
describes the integrable Lieb-Liniger model [5B] with en- 
ergies of different angular momentum states shifted with 
respect to each other by a Galilean transformation due to 
the rotational phase O. A finite barrier (b > 0) couples 
states with different total angular momentum, leading to 
the avoided level crossing seen in the inset in Fig. [T] At 
the precise position of the avoided crossing, an effective 
two-level system is realized with eigenstates being 50/50 
superpositions of the two angular momentum states. By 
adiabatically changing the applied phase from zero into 
the avoided crossing the superposition is created, while 
a rapid, non-adiabatic procedure leads to coherent oscil- 
lations between the two states. The oscillations with the 
period of AE/H, where AE is the level splitting at the 
avoided crossing [27], can be used to detect the presence 
of superposition states. The quantized total angular mo- 
mentum is measured by dropping the trap and observing 
the interference pattern develop a hole (K = Nh) or a 
peak (A" = 0) along the ring axis [19]. 

In order to observe coherent oscillations, the level split- 
ting AE should be larger than the rate of decoherence. 
We have calculated how AE scales with particle number 
in different regimes of interaction strength. While weak 
interactions lead to NOON states, the unfavorable scaling 
of AE and sensitivity to particle loss severely limit the 
attainable particle numbers. In the strongly-interacting 
regime, however, AE grows with particle number and the 
superposition states become robust against single parti- 
cle loss. 



III. NUMERICAL SIMULATION 

We have simulated the system for N = 5 particles 
by numerical diagonalization of the Hamiltonian (fijl. 
In order to obtain accurate results into the strongly- 
interacting regime, we have developed an effective Hamil- 
tonian for the one-dimensional Bose gas [28]. In Ap- 
pendix [A] we show that a rescaled interaction strength 
g = g(\ + g/go) 1 produces exact energies and eigen- 
state projections for two particles when the Hamiltonian 
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FIG. 2. (Color online) The level splitting AE between 
the ground and first excited state at Q — n is plotted over 
interaction strength g for 5 atoms with b/L = 0.008-Bo- 
The top horizontal axis shows the Lieb-Liniger parameter 
7 = g2it 2 /(EqLN). The dashes on the figure margins indi- 
cate analytic results for non-interacting (g = 0) and strongly 
interacting atoms (g — oo). Subplots show the distribution of 
the total angular momentum for (a) g — and (b) g — oo. 



is represented in a truncated occupation number basis 
from r angular momentum modes with g Q « rLE /2, 
where we have ignored small energy dependent terms. 
For larger particle numbers we use the same value of g 
and monitor the error by comparing with the exact Lieb- 
Liniger result. We still find a significant improvement 
compared to the unsealed Hamiltonian. We have used 
r — 20 modes for all calculations (except Fig. [3^,). In the 
Tonk-Girardeau regime, where the error is the largest, 
the relative error after rescaling is less than 3%, which 
is a factor of 8 times smaller than without the rescaling. 
In addition to the level splitting AE shown in Fig. [2] 
the simulations also provide insight into the nature and 
the robustness properties of the superposition states as 
shown in Fig. [3] 



IV. NON-INTERACTING ATOMS 

Before discussing the many-particle states in different 
interacting regimes we first consider the solutions for a 
single atom in the system (see Appendix [B]) . These are 
plane waves with angular momentum nh in the absence of 
the barrier. At Q = it the pairs of energy levels with ki = 
—nh and hi = (n+l)h are degenerate but this degeneracy 
is lifted for non-zero barrier strength b. Expressions for 
the energy levels e M and eigenstates are obtained using 
the transformation \&(x) = 4>(x)e lVtx / L , where <j>(x) 



x + L)e and employing the ansatz <f>(x) 



i2-nax/ L 



A p —i2-Kax/L 

For the case of $7 = tv, the discrete solutions a M with 
|U = 0,1,2,... correspond to the single particle energy 
levels e^ = a 2 Eo, and are the roots of 2wh 2 a ll /MLb = 



— tan(7ra M ) for odd \i and cfy = (/i + l)/2 for even /u. 
The i/th avoided crossing has a level splitting of E2 V -i — 
e 2(u-\) ~ V-^ f° r small barrier 6 -C 2vE L but reaches 
a constant value £2^-1 — £ 2(f-i) ^ (^ — 1/4) 2?o for large 
impenetrable barrier 6 > 2vEqL. 

For non-interacting atoms the energy gap A_E is the 
energy required to excite a single atom from the ground 
state across the first avoided level crossing (see dash on 
the far left of Fig. [2]) . Instead of a binary superposition, 
however, we find a binomial distribution involving many 
different momenta, as seen in Fig. [2^,. 



V. TONKS-GIRARDEAU REGIME 

The situation becomes more interesting when repulsive 
interactions remove the near degeneracies of all angular 
momentum eigenstates but the K — and the K = 
NH states. The interacting quantum system is generally 
very difficult to model and our numerical approach is 
limited to a small number of particles. One exception is 
the Tonks- Girardeau limit of strong interactions and low 
densities where 7 = gML/(h 2 N) ^> 1. Here, the bosons 
are strongly correlated as they cannot pass each other 
and undergo fermionization. The energy spectrum of the 
bosons in this limit is identical to that of non-interacting 
spinless fermions |29j but the single-particle momentum 
distribution is different and spreads oc l/-y/|fc| over an 
infinite range even at zero temperature |30j . Experiments 
have already validated the features of fermionization of 
ultra-cold atomic gases [20] [31] . 

A calculation of the energy gap in the Tonks- Girardeau 
limit proceeds as for non-interacting spinless fermions, 
using the single-particle results. For simplicity we as- 
sume an odd number of particles N. The ground state 
energy is found by summing the lowest N single-particle 
levels, while the first excited state energy is found by sub- 
stituting the highest level in the sum by the next higher 
level. The gap is given by the difference between these 
two energies AE — en — £jv— 1 an d is marked on the far 
right of Fig. [2] In the weak barrier regime determined by 
b < NE L, we find that 



AE = b/L + 0(b 2 ) 



(3) 



scales independently of particle number while the maxi- 
mum attainable gap energy is 



A£ max = (N/2 + l/A)E . 



(4) 



For a weak barrier and finite but strong interactions, 
the gap energy shows power law scaling with system 
size at constant density AE/(NEo) ~ N a with a = 
—4/7 + 0(7~ 2 ) according to to Luttinger liquid the- 
ory [35] [33], where NE$ is the energy of the first su- 
percurrent excitation. 

The distribution of the total angular momentum in the 
ground state is shown in the two subplots in Fig. [2] for 
(a) non-interacting atoms and (b) the Tonks- Girardeau 
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FIG. 3. (Color online) The quality of the superposition be- 
fore and after removal of an atom. The dashed line shows the 
superposition quality Q — 4P(0)P(Nh) for the ground state 
of the rotating system with parameters as in Fig. [2] The solid 
line shows the average superposition quality Q^ ' after the 
loss of an atom as described in the text. The crosses in in- 
set (a) show Q' -1 ' vs. atom number in the Tonks-Girardeau 
regime, and the solid line shows the result 1 — 1/N for non- 
interacting fermions for comparison. The horizontal bar vis- 
ible for N = 6 estimates the error in the numerical result 
from an insufficient number r of single particle modes by ex- 
trapolating from r — 14 (the reported result) to larger r by 
adding twice the difference from the r = 12 result. All other 
results use r — 20. Inset (b) shows the total angular momen- 
tum distribution after removal of an atom with k — lh in the 
Tonks-Girardeau regime. 



regime. The probability of the system having an an- 
gular momentum of K, the neutral atom equivalent of 
flux or current, is given by P{K) — Tr(Txp)- Here 
Tk = J2 n \Ki n ){Ki n \ ls the projector onto states with 
total angular momentum K, where the sum runs over all 
states with angular momentum K, and p = |^)(V'|/(V'|V') 
is the density operator. Subplot (b) clearly shows a 
superposition of two total angular momentum states, 
where, for N atoms, the total momentum difference of 
the two states is NH. We quantify the quality of a su- 
perposition between states with angular momentum Ki 
and A^ 2 by 



Q = 4P(A 1 )P( J K 2 ). 



(5) 



As desired, the superposition quality Q varies between 
the maximum value of 1 for an equal superposition and 
for angular momentum eigenstates. In Fig. [3] the dashed 
line shows Q for the ground state (with K\ = and 
K 2 = NH). We see that Q is very close to one beyond 
the minimum of AE at g « O.IEqL, where we also find 
that the distribution of total angular momentum does 
not change any more. The ground state is thus a binary 
superposition of total angular momentum for sufficiently 
large g. 



VI. NOON STATE REGIME 

The final feature of the AE curve in Fig. [2] that must 
be explained is the minimum around g ss O.IEqL. This 
is a regime where interactions are too weak to strongly 
correlate the atoms but strong enough to remove the de- 
generacies between states with intermediate angular mo- 
mentum. We can assume the ground state for Q = tt 
to take the NOON form \N,0) + \Q,N), where the two 
numbers in the ket denote the occupation number of the 
single-particle modes with angular momentum and h. 
For this to happen, the coupling to other angular momen- 
tum modes must be negligible. As shown in Appendix [C| 
this gives the conditions 



(6) 



bVN/L < gN/2L < E , 



where the mean interaction energy per particle is small 
compared to the kinetic energy per particle and large 
compared to the product of barrier energy and the square 
root of particle number. 

The energy splitting AE can now be calculated under 
the assumption that particles only access the two modes 
with k = and k = H. The Schrodinger equation is 
solved by eliminating states with intermediate angular 
momentum to yield (see Appendix [C[) , 



kN 



AE = 



N 



g N - 1 L(N-l)\ 



(7) 



We find that AE is small and decreases faster than expo- 
nentially with atom number N. This makes the NOON 
state experimentally difficult to prepare through tuning 
fi due to the long coherence times required [15] . This 
difficulty can be overcome by creating a superposition 
in the strongly-correlated regime first and then adiabat- 
ically reducing g. The time to reach the NOON state is 
then no longer limited by AE. The many-body dynamics 
involved in this process is left for future study. 



VII. PARTICLE LOSS 

Important limiting processes for the lifetime of ultra- 
cold atom experiments involve the loss of particles from 
the trap, e.g., due to collisions with high-energy atoms 
from the background gas. When the environment gains 
information about the state of the system this may lead 
to the collapse of the superposition. We consider the 
most detrimental case of information about the angu- 
lar momentum of the atom being gained by the envi- 
ronment. The resulting state after complete removal of 
a single atom with angular momentum kh from the N 
atom ground state |$) is modeled by 



We measure the robustness of the superposition state un- 
der particle loss by the averaged quality 

<5 hj =EgL"V/ jv "' (9) 



where 



Q 



4PI- 1 ! (~kh)P^ ((N - k)h) (10) 
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is averaged over the angular momentum kh of the lost 
particle weighted by the probability of finding an atom 
in mode k and n^ = (^JT '|aj,afc|^ _ )• Figured shows 
<5' _1 1 after single-particle loss alongside the Q for the 
ground state as a function of the interaction strength 
g. The regime around g sw O.IEqL features a high su- 
perposition quality but poor robustness against single- 
particle loss as expected for NOON states. In the Tonks- 
Girardeau regime, however, robustness increases dramat- 
ically and Q' -1 ' is still of the order of one! Inset (b) 
in Fig. [3] demonstrates that the binary superposition is 
maintained even after removal of an atom, while inset 
(a) shows how robustness increases with atom number. 
The increased robustness is consistent with the spread- 
out nature of the single-particle momentum distribution. 
Measuring a single-particle's momentum is not sufficient 
to determine the total angular momentum of the state. 



VIII. CONCLUSION 

We have presented a scheme for producing a robust 
binary superposition with strongly correlated atoms in 
the ground state of a rotating system. A barrier couples 
two states that differ by angular momentum Nh where 
N is the number of atoms and the maximally achievable 
level splitting is proportional to N. Due to this favorable 
scaling, mesoscopic superpositions involving hundreds or 
thousands of atoms become feasible. Loading, e.g., 100 
atoms of 7 Li into a ring trap with radius 50/mi leads 
to a mean particle spacing of 3.1/xm. An energy gap of 
AE rj 25Eq rj 45/iHz (half the limiting value) could be 
realized in the Tonks-Girardeau regime with a barrier 
rotating at angular speed u = 0.29 x 2ttKz. Such states 
are much less fragile than the related NOON states and 
we expect a robustness of Q' -1 ' > 1 — 1/JV = 0.99, based 
on extrapolating from Fig. [3k. 
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Appendix A: Effective Hamiltonian for numerical 
simulations: 



For the numerical calculations we perform an exact di- 
agonalization of the Hamiltonian in a truncated Hilbert 
space. In order to improve the accuracy of the trun- 
cated calculation we have developed an effective Hamil- 
tonian that involves rescaling the interaction constant 
as detailed in the following. We start with the second 
quantized form of the Hamiltonian H = Hk + Hb + Hi, 
where Hk, Hb, and Hi are the Hamiltonians describ- 
ing the kinetic energy of the atoms, the barrier, and the 
interactions between the atoms, respectively. These are 
given by 
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Hr 



Ht = 



dx b 8(x) &(x)4>(x) 



dx I ¥(x)&(x)$(x)i>{x), 



(Al) 



where &'(x) and ^(x) are the Schrodinger field opera- 
tors with standard bosonic commutation relations. We 
transform into a truncated angular momentum basis with 
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1 
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(A2) 



where d k and a^ create and destroy an atom with angular 
momentum kh, respectively. We introduce the effective 
Hamiltonian in the truncated basis by H — Hk+Hb+Hi 
with 
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H, 
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E 4 I «L 6fc i-? a *2+9' ( A3 ) 



fci,*3,g=-r/2+l 



which becomes formally exact and identical to H of 
Eq. (Al) for r — > oo and g = g. Here, r is the num- 



ber of angular momentum modes used in the simulation, 
which we choose to be even. Choosing finite r effectively 
truncates Hilbert space and constitutes an approximation 
that is expected to converge towards the exact result for 
large r. We have used up to r = 36 for testing conver- 
gence. The results presented in Figs. [51 and [3] have been 
obtained with r = 20 (except for Fig. km). While at this 
level the ground and excited state energies have already 
converged to the level of machine precision for small g 
up to the NOON-state regime (g sw 0.1), we find slow 
convergence with increasing r in the strongly-interacting 
Tonks- Girardeau regime, where we can compare with ex- 
act results from the Bose-Fermi mapping [29 . 



Using an effective Hamiltonian for our numerical cal- 
culations is a way to improve the accuracy significantly. 
Formally, the effective Hamiltonian is introduced by a 
transformation of the full Hamiltonian onto the truncated 
Hilbert space that preserves a subset of the exact eigen- 
values [28] . Here, we use a convenient approximate form 
that involves solely a rescaling of the interaction constant 
g in Eq. ( A3 1 , found by considering the simple system of 
two interacting particles at b — and fl = 0, which we 
solve analytically. 

We look for the ground state of H for two par- 
ticles and write the general wave function as \ip) — 



fcl, 



Cfci.fe^LaL 



vac). Substituting this into the 



Schrodinger equation and making use of bosonic symme- 



try C fcl>fc2 
equations 



= C, 



fco./i'l 



we obtain the set of simultaneous 



{E-E Q {kj + k 2 2 ))C klM = ^ J2 C ' 



-q,k 2 



9 (A4) 



All terms in the Hamiltonian conserve total angular mo- 
mentum and we can expect the ground state to have zero 
angular momentum where k\ = —k 2 . It follows that the 
right hand side of Eq. ( A4 ) is independent of k\ . Using 

oc (E-2E q 2 y 1 and 



this it is easy to show that C- qq 

L °° L 
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q— — oo 

where L/g = Eg=_ r / 2 - 

would obtain from the truncated Hamiltonian H 

Eq. (Al). The term -L/g 



n 



is the result that we 

of 
accounts for the sum over 
the remaining angular momentum modes and thus de- 
fines <7o- From rearranging Eq. ( A5 1 we find that both the 
energy E as well as the coefficients C- qA obtained from 
diagonalizing the effective Hamiltonian H agree with the 
exact results if the interaction strength is rescaled to 



= .9/(l + 5/5o), 



where 50 is given by 
L _ 2 
go E r 



2E- E 
6E 2 r 3 



+ 0(r- 5 ). 



(A6) 



(A7) 



The terms beyond the first on the right hand side explic- 
itly depend on the energy E of the solution. Due to their 
scaling with the number of single-particle modes r, these 
terms can be neglected when r 2 Eo ^S> E — Eq. For more 
than two particles the formally exact effective Hamilto- 
nian contains three and more particle interaction terms, 
which we will ignore here |28j . For the calculations in 
this work, we approximate the effective Hamiltonian by 
H of Eq. |A3l with g = rLE /2. 



Appendix B: Single particle spectrum: 

The Hamiltonian describing the system of a ID loop 
with a rotating barrier is given in Eq. ([I]). For a single 



atom analytic solutions can be found. We write the wave 
function as ty(x) — 4>(x)e inx / L and substitute this into 
the Schrodinger equation to obtain 

-^-^4>{x) + b5(x)4>(x)=e<f>(x), (Bl) 

where the boundary conditions require <f>{x) = <f>(x + 
L)e lQj . The first derivative of the wave function at the 
Dirac delta barrier is discontinuous and is found by inte- 
grating the Schrodinger equation over the barrier: 



dx 



,,'H 



x=+0 



dx 



2Mb 



0(0). 



(B2) 



x=L-Q 

With the ansatz cf>(x) = e l27rax / L + Ae~ i2 ™ x / L and us- 
ing the boundary condition for the wave function we find 



.4: 



i27ra 



S and S = sin(vra + fi/2)/ sin(7ra-£>/2). Sub- 



stituting this into Eq. ( B2 ) yields 
cot(7ra 



£5/2) 



-cot(7ra + £>/2), (B3) 
The discrete solutions a,. 



47Tft a 

MLb 
which has to be solved for a. 
with fi = 0,1,2,... correspond to the single-particle en- 
ergy levels e M = a^Eo. For the case of O = x, Eq. ( B3 ) is 
simplified to 2Trti 2 a ll /MLb 
ttp = (// + l)/2 for even /i. 



•tan(Tra^) for odd fi and 



Appendix C: Conditions for existence and energy 
gap for NOON state: 

The values of the experimental parameters needed to 
create a NOON state in this system are determined by 
three factors [16]. Firstly we want the states |0, N, 0) 
and |0, 0, N) to be degenerate (here we describe the state 
of the system in the occupation number basis where the 
three numbers in the ket represent the —h,0 and h an- 
gular momentum modes, respectively). This is achieved 
when n = 7r, which is the condition we are considering 
here. 

Secondly, the two states must be, at most, weakly 
coupled to and energetically well separated from other 
states, otherwise the ground state will be a superpo- 
sition of many states. This is achieved by requiring a 
minimum interaction strength to energetically separates 
the states |0,iV — n,n), where n = 1,...,N — 1 from 
|0, iV", 0) and |0, 0, N). At the same time, the interaction 
should be weak enough to not couple strongly to other 
states through interactions. We can estimate these cri- 
teria by considering two simple two level systems. First 
we consider the states |0, iV, 0) and \0,N — 1,1), where 
\0,N — 1,1) is the state that is coupled strongest to 
|0, N, 0) through the barrier term in the Hamiltonian. 
The Hamiltonian for this two state system is, 

# e = f(iV-l)|0,iV-l,l)<0,iV-l,l| 



A(|0,iV,0)(0,A-l,l| 

+|0,JV-l,l)(0,iV,0|) (CI) 



where we have ignored a constant energy term. Starting 
with the ansatz |\E>) = a \0, N,Q) +Oi|0, N- 1, 1) we find 




2bVN 



2b\fN 



1. 



(C2) 



We must also consider the relative amplitudes of the 
states |0,A,0) and |1,JV-2,1), where 1 1,7V -2,1) is the 
state that is coupled strongest to |0, N, 0) through the 
interaction term in the Hamiltonian. The Hamiltonian 
for this two state system is, 

H e = \2E + -(2N-3)] |1,JV-2,1)(1,JV-2,1| 

+ ^N(N-l)(\0,N,0)(l,N-2,l\ 

+ |1,A-2,1)(0,A,0|) (C3) 

Again we have ignored a constant energy term. This time 
we start with the ansatz |*) = a |0, N, 0)+5i|l, N — 2, 1) 
and find, 



2E L + g(2N - 3) 
g^/N(N-lT~ 



{2E L + g{2N - 3) 



(C4) 



To create a NOON state we require |ao/ai|, |ao/ai| >• 1, 
which is achieved when, 



(C5) 



'N qN 

L «-2L <<E °- 



We see that the mean interaction energy per particle 
needs to be much smaller than the kinetic energy per par- 
ticle and much larger than the barrier energy times the 
square root of the number of particles. Eq is inversely 
proportional to the square of the circumference of the 
loop, L 2 , and the mass of the atoms, M. Therefore, it 



is easier to fulfill condition (C5) with lighter atoms and 



smaller rings. As the number of particles are increased it 
will be experimentally unattainable to reach the small in- 
teraction strength and the barrier height needed to make 
a NOON state. 

Finally, in order to obtain the superposition, the states 
|0, N, 0) and |0, 0, N) must be coupled. By calculating the 
coupling strength we will also be able to calculate the en- 
ergy level splitting between the ground and first excited 
states. There is no direct first order coupling between 
the two states, but they do couple through intermediate 
states. We have already argued that the population in 
states other than |0, TV, 0) and |0, 0,iV) is small. In this 
regime and at £2 = n it is a good approximation to con- 
sider just the and 1 angular momentum modes. States 
that have atoms with other angular momentum modes 
have a larger kinetic energy associated with them. This 
makes the states energetically unfavorable and only pro- 
vide a small addition to the coupling strength between 
|0,JV,0) and |0,0, AT). 
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FIG. 4. (Color online) NOON state energy level splitting as 
a function of atom number. The line shows the energy level 
splitting calculated using Eq. fty, where b/L — 0.008i?o and 
g/L = 4nby^N/L(N - 1). In the regime where the NOON 
state is formed it is adequate to describe the system using 
just the and ft angular momentum modes. The results from 
the numerical simulation are shown by the red crosses. The 
energy level splitting decreases rapidly as the number of par- 
ticles are increased. We see the numerical result break down 
for TV > 11 due to limited numerical accuracy. The analytic 
result is still valid beyond this point. 



The Schrodinger equation for this system can be writ- 
ten as a set of simultaneous equations, 

Xa n = t n a n + Vn t n-\a n -i + V nt n+\0^n+li (C6) 

Here a n is the coefficient of state \0,N — n,n), A are 
the eigenenergies that solve the system, t n — (0, N — 
n, n\H\0, N — n,n) — gn{N — n)/L + const., V n , n +\ = 
(Q,N-n,n\H\0,N-n-l,n+l) = b y/(N -n)(n + Tj/L 
and H is the Hamiltonian. We can systematically elim- 
inate coefficients of intermediate states leaving just arj 
and ajv, 



a N = 



(A-t )(A-*i)...(A-tjv-i) 



VoiVw-Vn- 



A 



N 



1,N 



a (C7) 



To prove by induction that this is the general form we 
add another atom to the system and eliminate a^r using 
the additional equation Xa^ = t^a-N + Vjv\jv-iGjv-i + 
Vn,n+io,n+i- This leaves just arj and qn+i and gives 





r(A-t )(A-*i)...(A-tjv) ! 1 




&N+1 


^01^L2---VAf,JV+l 


a 0i 


where 


A A {X ' tN) 

An+i — An 

Vnm+i 




Vn,n-i (A - t )(A - ti)...(A - t N _ 2 ) 


Vn,N+1 Vq{Vi2 — Vn-2,N-1 


, Vn.n-i 
A N -! 

VN.N+l 

and A x = 0. 


We are 


now left with two simultaneous equations 



(C9) 



Aa = t(X)a + V(X)cin, 
Xa N — t(X)a N + V(X)a , 



(C10) 



where 



V m V 12 ...Vi 



V{X) (X-h)(X-t 2 ). 
t(X) = t -A N V(X). 



N-1,N 



.(A-ijv-t)' 



(Gil) 



When the amplitude of states other than |0,iV, 0) and 
|0, 0, N) are small, Eq. (C10) gives the ground and first 
excited state energies, because |0, N, 0) and |0, 0, N) have 
the lowest energies of the uncoupled system. If we assume 
i(A ) w t(Ai) w t(t ) and V(X ) « V(Ai) « V(t ) then 
we have a simple two level system and the energy level 
splitting is given by AE — Ai — Ao = 2V. We obtain 
Eq. tf?\. For large N this expression is dominated by 
the inverse factorial and thus becomes exceedingly small. 
We have verified the assumptions by expanding around t 
and find it is valid under the conditions given in Eq. ( C5 1. 
We compare Eq. Q with a full numerical solution of the 
two mode model of Eq. ( C6 ) in Fig. BJ 
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